Automatic implementation of material laws: Jacobian calculation in a finite element code with TAPENADE

نویسندگان

  • Florian Zwicke
  • Philipp Knechtges
  • Marek Behr
  • Stefanie Elgeti
چکیده

In an effort to increase the versatility of finite element codes, we explore the possibility of automatically creating the Jacobian matrix necessary for the gradient-based solution of nonlinear systems of equations. Particularly, we aim to assess the feasibility of employing the automatic differentiation tool TAPENADE for this purpose on a large Fortran codebase that is the result of many years of continuous development. As a starting point we will describe the special structure of finite element codes and the implications that this code design carries for an efficient calculation of the Jacobian matrix. We will also propose a first approach towards improving the efficiency of such a method. Finally, we will present a functioning method for the automatic implementation of the Jacobian calculation in a finite element software, but will also point out important shortcomings that will have to be addressed in the future.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Fast Automatic Differentiation Jacobians by Compact LU Factorization

For a vector function coded without branches or loops, a code for the Jacobian is generated by interpreting Griewank and Reese’s vertex elimination as Gaussian elimination and implementing this as compact LU factorization. Tests on several platforms show such a code is typically 4 to 20 times faster than that produced by tools such as Adifor, Tamc, or Tapenade, on average significantly faster t...

متن کامل

Fast Ad Jacobians by Compact Lu Factorization

For a vector function, coded without branches or loops, code for the Jacobian is generated by interpreting Griewank and Reese’s Vertex Elimination as Gaussian elimination, and implementing this as compact LU factorization. Tests on several platforms show such code is typically 4 to 20 times faster than that produced by tools such as Adifor, Tamc or Tapenade, on average significantly faster than...

متن کامل

Buckling Analysis of Rectangular Functionally Graded Plates with an Elliptic Hole Under Thermal Loads

This paper presents thermal buckling analysis of rectangular functionally graded plates (FG plates) with an eccentrically located elliptic cutout. The plate governing equations derived by the first order shear deformation theory (FSDT) and finite element formulation is developed to analyze the plate behavior subjected to a uniform temperature rise across plate thickness. It is assumed that the ...

متن کامل

FINGER: A Symbolic System for Automatic Generation of Numerical Programs in Finite Element Analysis

FINGER iS a Lisp-based system to derive formulas needed in finite element analysis, and to generate FORTRAN code from these formulas. The generated programs can be used with existing, FORTRAN-based finite element analysis packages. This approach aims to replace tedious hand computations that are time consuming and error prone. The design and implementation of FINGER are presented. Techniques fo...

متن کامل

Forward-Mode Automatic Differentiation in Julia

We present ForwardDiff, a Julia package for forward-mode automatic differentiation (AD) featuring performance competitive with low-level languages like C++. Unlike recently developed AD tools in other popular high-level languages such as Python and MATLAB, ForwardDiff takes advantage of just-in-time (JIT) compilation to transparently recompile AD-unaware user code, enabling efficient support fo...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:
  • Computers & Mathematics with Applications

دوره 72  شماره 

صفحات  -

تاریخ انتشار 2016